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Abstract 

We resum the leading ultrasoft logs of the singlet and octet static QCD potentials 
within potential NRQCD. We then obtain the complete three-loop renormalization 
group improvement of the singlet QCD static potential. The discrepancies between 
the perturbative evaluation and the lattice results at short distances are slightly 
reduced. 
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Recently, there has been quite a lot of effort in identifying the large logs that arise from 
the different scales involved in heavy quark systems near threshold [1-7]. The proper renor- 
malization group (RG) rcsummation of these logs is a non-trivial issue. It has recently been 
addressed within the so-called vNRQCD approach [8] for the 0(l/m^) potentials [9], and 
for the 0{l/m) potential and the production current [10]. Since the proper resummation of 
the logs may be important for the physics of the top quark production near threshold, as 
well as for heavy quarkonium systems, and the existing calculations (see [2]) arc done in a 
formalism closer to potential NRQCD (pNRQCD) [11,4], it would be desirable to know how 
to RG-improve within the latter formalism. In this note we provide the first step towards 
this goal by showing how the static potentials of QCD can be RG-improved within that 
framework. Moreover, the RG-improved static potentials obtained in this paper represent a 
new result by themselves. They provide the complete three-loop RG evolution of the static 
potentials within an expansion in as (in order to be so, besides our calculation, one needs 
to know the static potentials at two loops [12,13]). Our results may also be relevant in or- 
der to understand the discrepancies, at relatively short distances, between the perturbative 
evaluation and the lattice results. 

Since we are only interested in the static potentials, we only need to consider the static 
limit of NRQCD [14] and pNRQCD, i.e. we only need to work at leading order in 1/ni. 
The matching between NRQCD and pNRQCD in the static limit, in the situation where 
Aqcd <^ 1/r (the limit we will consider in this paper), has been worked out in detail in ref. 
[4]. 

The pNRQCD lagrangian at leading order in 1/m and next-to- leading order in the multipole 
expansion reads 



>CpNRQCD = Tr jst (ido -V,)S + Qt (iDo - K) o| 



+gVA{r)Tr {O^r • E S + S^r • E O} + c/^^^^Tr {O^r • E O + O^Or • e} 
1 



F''F^'"'^0{r'). (1) 

All the gauge fields in Eq. (1) are evaluated in R and in particular F^^°- = F'*''"(R, t) and 
iDoO = idoO - g[Ao{R, t), O] and 

The potentials Vj, i = s,o, A, B are to be regarded as matching coefficients, which depend on 
the scale u separating soft gluons from ultrasoft ones. In the static limit, we understand by 
soft energies the ones of 0{l/r) and by ultrasoft energies the ones of 0{as/r). Notice that the 
hard scale, m, plays no role in the static limit. The only assumption made so far concerning 
the size of r is that 1/r ^ Aqcd- However, in order to use perturbative RG techniques, we 
shall also assume from now on that we are working at scales u such that as{i') <^ 1. 
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Formally, we can write Eq. (1) as an expansion in r in the following way: 

oo 



n=-l 



where the above fields and parameters should be understood as bare and the renormalization 
group equations of the renormalized matching coefficients read 

u^\ = B,{X). (3) 

If A_i = 0, there are no relevant operators (superrenormalizable terms) in the Lagrangian 
and the RG equations have a triangular structure (the standard structure one can see, for 
instance, in HQET [15]; see Ref. [16] for a review): 

u-^Xo = So(Ao) (4) 
i/-^Ai = Si(Ao)Ai (5) 
i^^A2 = -62,(2,1) (Ao)A2 + -82,(1,2) (Ao)A^ (6) 



(7) 



where the different B's can be power-expanded in Aq (Aq corresponds to the marginal op- 
erators (renormalizable interactions)). If A_i ^ 0, however, there are relevant operators 
(superrenormalizable terms) in the Lagrangian and the RG equations lose the triangular 
structure. Still, if A_i <^ 1, a perturbative calculation of the renormalization group equa- 
tions can be achieved as a double expansion in A_i and Aq. The RG equations now have the 
following structure: 

= S_i(Ao)A_i + S(_i,2)(Ao)A^Ai + 3)(Ao)A^A? + 3^(Ao)AliA2 + O(Aii) 

z/^Ao = So(Ao) + So,i(Ao)A_iAi + O(A^) (9) 

u^X, = B^{Xo)X, + Si,i(Ao)A_iA? + 0{X\) (10) 
and so on (for a similar discussion in the context of scalar A0"-like theories see [17]). 
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At short distances, the static hmit of pNRQCD hves in the second situation. Specifically, 
we have A_i = {ay^,av<,}, that fulfils A_i ^ 1; Aq = and Ai = {Va, Vb}- Therefore, we 
can calculate the anomalous dimensions order by order in as- In addition, we also have an 
expansion in A_i that corresponds to working order by order in the multipole expansion. 

The specific form of the pNRQCD lagrangian severely constrains the above general structure. 
From ref. [4]Q we obtain at leading, non-vanishing, order 



u-^as = asP{as) 
dv 

y^VA = Q + 0{\_,\oAl) (11) 
i.Ay^ = + 0(A_iAo,A2). 



It is the aim of this work to solve Eq. (11), and hence to provide the leading log (LL) ultrasoft 
RG improvement of the pNRQCD lagrangian in the static limit. 

The last two equations in Eq. (11) are decoupled from the rest and are equal to zero at the 
order we are working. Therefore, we immediately obtainQ 

VA{y)=VA{r) and ^/^(z/) = ^^^(r). (12) 



The first two equations of Eq. (11) are equivalent to 



J^-i-{(^Vo-(^Vs) = 0, (13) 



du 



where 



7o. = ^. (14) 

^ The result of the last two equations corrects the previous evaluation of Ref. [4] (we thank A. Vairo 
for confirming these results). In order to cross-check this result we have redone the computation 
both in the Coulomb gauge and in the background field gauge. 

^ Throughout the paper, we use as{r) as a shorthand for as{i^ = 1/?") and similarly for any other 
running object at the scale v = 1/r. 
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The second equation tells us that the combination ay^ — c^Vs is RG-invariant at this order. 

The first equation can now be easily solved (note that in order to solve these equations with 
the demanded accuracy it was necessary to know the RG solution of V^, which happened to 
be trivial). We obtain (where /3o = llC^/3 — 4Tpn//3, Tp = 1/2 and n/ is the number of 
light quarks): 



X 



1 - 



V 



1 _ ±Za 



^ ((^ - Cf) avXr) + CfavXr)) VI{t) log (f^) j 
Ca 



37o 



I 



X 



(15) 



which completes our results. These are exact up to O \ y?_Xj) (^A^]^(r) (Ao(r) logrz/)"j , 

A_i(r)Ao(r) (^A^ ^(r) (Ao(r) log rz/)"^ ^ where n = 0, 1, 2... and m = 1, 2... . At this point we 
would like to stress the simplicity of the calculation, which follows to a large extent from the 
formalism used. 



In the above results the different matching coefficients can be considered to be independent 
(though assuming A_i, Aq ^ 1 in order for our perturbative evaluation to make sense). If we 
want to perform a strict expansion in Og, one has to use the perturbative relation between 
A_i(r), Ai(r) and a^ij'): A_i(r) = A\ajyr) + A2Q;s(r) + ... and = 1 + Bia^r) + .... Since 

ay^ and ay^ are known at two-loop accuracy [12,13], Eq. (15) is known with the following 
accuracy 

OiyXv)^Oi%]^'''"'\r) + ^fAr) log + 0(a,(r)^+" log" r^) (16) 

avXy)-olil^''^'\r) + ^al{r) log + 0(«,(r)^+" logV^.), 

where n = 0, 1, 2, ... . 

The above equation represents the complete three-loop RG improvement of the static po- 
tentials. When the running Q;s(r) is substituted by its three-loop expression above and ex- 
panded, all the log rz/ in the three-loop static potentials are obtained, including the infrared 
logs of [18,3]. More explicitly, if we write A_i = A_i(r, aXy))^ Eq- (16) provides the correct 
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Fig. 1. We plot ay^ versus r (GeV) for a range of energies relevant to a comparison with lattice 
simulations and also to the study of the T(15) system (nj- = 4, Aq^'^j^ = 0.28 GeV). The 

solid line corresponds to a^^^""^^"*, the dotted line to oty^'^'^^^^ plus the single leading-log correction 
computed in Ref. [3], the dot-dashed line to Eq. (16) and the dashed line to Eq. (15). We have fixed 
V = as{r)/r. 

Q;g+"'(z/) log" rz/, ag"''"(i^) log" rz/, ^^"''"(z/) log" rz/ terms at any power of n. Note that such 
accuracy is not achieved by only replacing as{r) by its three-loop running expression in the 
two-loop calculation of the potential, as it is sometimes done in the literature [19] (see also 
[20]). 

The Lagrangian (1) is suitable for the study of static systems. In particular, we may be inter- 
ested in computing the energy of two infinitely heavy (static) sources separated at a distance 
r in a singlet configuration. This is, indeed, the object computed in lattice simulations and 
usually referred as the static potential. Within our formalism, this energy corresponds to 
the pole of the singlet- singlet correlator. In the energy, the scale dependence u oi Vg gets 
cancelled by ultrasoft contributions of 0{as{r)/r). Therefore, by lowering the scale v in Vs to 
values of 0(as(r)/r), all the large logs arising from the ratio between the soft and ultrasoft 
scale are resummed and the energy becomes, at the accuracy we are working, equal to Vg 
with V = as{r)/r (up to finite terms). Now we would like to give some numerical estimates of 
our results (we fix z/ = as(r)/r in what follows). In Fig. 1, we plot ay^ versus r (GeV), within 
different approximations. If we compare the last term in Eq. (16) with the single leading-log 
result obtained in [3], we can see a sizeable correction of the same order of magnitude. If 
for definiteness we choose r^^ = 2 GeV, Eq. (16) gives a correction of around 3%, whereas 
the single leading-log term gives a correction of around 1%. This shows the necessity to 
RG-improve the single leading logs. In any case, as we can see from Fig. 1 and the above 
results, the correction is rather tiny whith respect to the total value of avg- Equation (15), 
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however, produces a much more sizeable correction, as we can see in Fig. 1. If we again take 
— 2 GeV the correction is of around 11%. 

In all the cases the new results lower the value of ay^, reducing the existing discrepancies 
between the perturbative results and the lattice simulations [21]. This is specially so when 
using Eq. (15). Although the good behaviour of ay^ using Eq. (15) is certainly appealing, 
wc should bear in mind that around the 1 GeV region A_i(r) ~ 1 and hence perturbation 
theory becomes unreliable. This is so specially in the perturbative relation between a.v^{r) 
and as{r) and it is usually associated to renormalon effects [22]. This point should be better 
understood before a detailed study about the impact of the RG improvement can be done. 

Since for «yo(r) only the one-loop expression is available, we refrain from performing a similar 
study for it (note that in the numerical analysis above we have used the one-loop expression 
of av„{r) in Eq. (15)). 

One could also consider non-perturbative effects in the static case when 1/r ^ Aqcd (see 
[4] for a more general discussion). Here we only mention that when as{r)/r 3> Aqcd the 
leading non-perturbative contribution is proportional to r^{asFF) /as{r) [24], whereas for 
Ois{f)/f ~ Aqcd one has a non-perturbative contribution proportional to times a non-local 
correlator of two chromoelectric fields with multiple potential insertions of 0{as{r)/r). 

Our results are also relevant to actual physical systems composed of a quark and an an- 
tiquark with very large but finite mass. However, in order to achieve, for instance, an 
0(TOQ;g"'""log"' QJs) accuracy in the binding energy of those systems, the RG improvement 
of the 1/m and 1/m^ potentials is also necessary. Since local field redefinitions allow to 
reshuffle contributions from a given order in 1/m to another [23], a meaningful outcome can 
only be obtained when all the contributions from the various potentials are taken into ac- 
count. Hence, we restrict ourselves here to present a few numerical estimates on the impact 
of the RG-improved static potentials for some typical scales one can find in the T(15') and t-t 
systems. In particular, for the T(15'), for which typical soft scales are of 0(2 GeV), the dis- 
cussion of Fig. 1 applies. For t-i systems, in order to give an estimate, we take = 20 GeV 
{uf = 5, Aq^q'""^*^ = 0.2). We obtain that the corrections are of order 0.4%, 0.5%, 1%, if we 
consider the single leading-log result, Eq. (16) and Eq. (15), respectively. 

Before closing, let us mention that in the vNRQCD approach it appears to be crucial that 
ultrasoft and soft gluons run from a scale mv"^ and mv, respectively, to the scale m. In the 
static system, the scale m does not exist, and hence it is not clear to us how one should 
proceed in the above-mentioned approach. Notice also that the difficulties pointed out in [25] 
concerning a nai've RG improvement in two stages do not apply to the static case, since here 
only one stage is involved, namely the RG improvement between NRQCD and pNRQCD. 

Acknowledgements 

We thank M. Beneke, J. Kuti, A. Manohar, I. Stewart and A. Vairo for stimulating dis- 
cussions. We also thank G. Bali for sharing his lattice data with us. A. P. acknowledges the 



6 



TMR contract No. ERBFMBICT983405 and J.S. the financial support from the CICYT 
(Spain), contract AEN98-0431, the CIRIT (Catalonia), contract 1998SGR 00026, and the 
grant BGP-08 (Catalonia). J. S. thanks A. Manohar and the HEP group at UCSD for their 
warm hospitality while this work was carried out. 



References 

[1] M. Beneke, A. Signer and V.A. Smirnov, Phys. Rev. Lett. 80 (1998) 2535; A. Czarnecki and 
K. Melnikov, Phys. Rev. Lett. 80 (1998) 2531. 

[2] A.M. Hoang, M. Beneke, K. Melnikov, T. Nagano, A. Ota, A.A. Penin, A.A. Pivovarov, A. 
Signer, V.A. Smirnov, Y. Sumino, T. Teubner, O. Yakovlev and A. Yelkhovsky, Eur. Phys. J. 
direct C 3 (2000) 1. 

[3] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D 60 (1999) 091502. 
[4] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566 (2000) 275. 
[5] B. A. Kniehl and A. A. Penin, Nucl. Phys. B 563 (1999) 200. 
[6] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Lett. B 470 (1999) 215. 



[7] B. A. Kniehl and A. A. Penin, |hep-ph/9911414 . 



[8] M. E. Luke, A. V. Manohar and I. Z. Rothstein, Phys. Rev. D 61 (2000) 074025. 
[9] A. V. Manohar and I. W. Stewart, |hep-ph/991222^ ; |hep-ph/000303^ . 



[10] A. V. Manohar and I. W. Stewart, |hep-ph/000310^ . 

[11] A. Pineda and J. Soto, Nucl. Phys. Proc. Suppl. 64 (1998) 428. 

[12] Y. Schroder, Phys. Lett. B 447 (1999) 321; M. Peter, Nucl. Phys. B 501 (1997) 471. 
[13] Y. Schroder, private communication. 

[14] W. E. Caswell and G. P. Lepage, Phys. Lett. B 167 (1986) 437; G.T. Bodwin, E. Braaten and 
G.P. Lepage, Phys. Rev. D 51 (1995) 1125, Erratum ibid. D 55 (1997) 5853. 

[15] S. Balk, J.G. Korner and D. Pirjol, Nucl. Phys. B 428 (1994) 499; C. Bauer and A. V. Manohar, 
Phys. Rev. D 57 (1998) 337; C. Balzereit, Phys. Rev. D 59 (1999) 034006. 

[16] M. Neubert, Phys. Rep. 245 (1994) 259. 

[17] M. Atance and J. L. Cortes, Phys. Rev. D 56 (1997) 3611. 

[18] T. Appelquist, M. Dine and I. J. Muzinich, Phys. Rev. D 17 (1978) 2074. 

[19] M. Peter, Phys. Rev. Lett. 78 (1997) 602. 



[20] V. V. Kiselev, A. E. Kovalsky and A. I. Onishchenko, |hep-ph/000502C 



7 



[21] G. S. Bali and K. Schilling, Phys. Rev. D47 (1993) 661; G. S. Bali and K. Schilling, Int. J. 
Mod. Phys. C4 (1993) 1167; G. S. Bah, K. Schilling and A. Wachter, Phys. Rev. D56 (1997) 
2566; K. J. Juge, J. Kuti and C.J. Morningstar, |hep-lat/9911007 . 

[22] M. Jezabek, J.H. Kuhn, M. Peter, Y. Sumino and T. Teubner, Phys. Rev. D 58 (1998), 014006; 
M. Jezabek, M. Peter and Y. Sumino, Phys. Lett. B 428 (1998) 352. 

[23] N. BrambiUa, A.Pineda, J. Soto and A. Vairo, |hep-ph/0002250[ 

[24] C.A. Flory, Phys. Lett. B 113 (1982) 263; R.A. Bertlmann and J.S. Bell, Nucl. Phys. B 227 
(1983) 435. 



[25] A. V. Manohar, J. Soto and I. W. Stewart, |hep-ph/000609e. 



8 



